Magnetic-field influenced non-equilibrium transport through a quantum ring 
with correlated electrons in a photon cavity 
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We investigate magnetic-field influenced time-dependent transport of Coulomb interacting elec- 
trons through a two-dimensional quantum ring in an electromagnetic cavity under non-equilibrium 
conditions described by a time-convolutionless non-Markovian master equation formalism. We take 
into account the full electromagnetic interaction of electrons and cavity photons without resorting to 
the rotating wave approximation or reduction to two levels. A bias voltage is applied to semi-infinite 
leads along the x-axis, which are connected to the quantum ring. The magnetic field is tunable to 
manipulate the time-dependent electron transport coupled to a photon field with either x- or y- 
polarization. We find that the lead-system-lead current is strongly suppressed by the y-polarized 
photon field at magnetic field with two flux quanta due to a degeneracy of the many-body energy 
spectrum of the mostly occupied states. Furthermore, the current can be significantly enhanced by 
the y-polarized field at magnetic field with half integer flux quanta. 
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I. INTRODUCTION 

Quantum interference phenomena are essential when 
developing quantum devices. Quantum confined geome- 
tries conceived for such studies may consist of which- 
path interferometers, 1,2 coupled quantum wires, 3,4 side- 
coupled quantum dots, 5,6 or quantum rings. 7,8 These 
coupled quantum systems have captured interest due 
to their potential applications in electronic spectroscopy 
tools 9 and quantum information processing. 10 Further- 
more, the magnetic flux through the ring system can 
drive persistent currents 11 and lead to the topologi- 
cal quantum interference phenomenon known as the 
Aharonov-Bohm (AB) effect. 12 16 Both, the persistent 
current and ring conductance show characteristic oscilla- 
tions with period of one flux quantum, <£ = hc/e. Vary- 
ing either the magnetic field or the electrostatic confining 
potentials allows the quantum interference to be tuned. 17 

There has been considerable interest in the study of 
electronic transport through a quantum system in a 
strong system-lead coupling regime driven by periodic 
time-dependent potentials, 18 21 longitudinally polarized 
fields, 22 24 or transversely polarized fields. 25,26 On the 
other hand, quantum transport driven by a transient 
time-dependent potential enables development of switch- 
able quantum devices, in which the interplay of the elec- 
tronic system with external perturbation plays an impor- 
tant role. 27 30 These systems are usually operated in the 
weak system-lead coupling regime and described within 
the wide-band or the Markovian approximation. 31 33 
Within this approximation, the energy dependence of the 
electron tunneling rate or the memory effect in the sys- 
tem are neglected by assuming that the correlation time 
of the electrons in the leads is much shorter than the 
typical response time of the central system. However, 
the transient transport is intrinsically linked to the co- 
herence and relaxation dynamics and cannot generally be 



described in the Markovian approximation. The energy- 
dependent spectral density in the leads has to be included 
for accurate numerical calculation. 

In order to explicitly explore the transport dynamics 
with transient system-lead coupling and electron-photon 
coupling, a non-Markovian density-matrix formalism in- 
volving the energy-dependent coupled elements should 
be considered based on the generalized master equa- 
tion (GME). 34 37 How to appropriately describe the car- 
rier dynamics under non-equilibrium conditions with re- 
alistic device geometries is a challenging problem. 38,39 
More recently, manipulation of electron-photon coupled 
quantum systems embedded in an electromagnetic cav- 
ity has become one of the most promising applications 
in quantum information processing devices. Utilizing 
the giant dipole moments of inter-subband transitions in 
quantum wells 40,41 enables researchers to reach the ul- 
trastrong electron-photon coupling regime. 42 44 In this 
regime, the dynamical electron-photon coupling mech- 
anism has to be explored beyond the wide-band and 
rotating-wave approximations. 45 47 Nevertheless, time- 
dependent transport of Coulomb interacting electrons 
through a topologically nontrivial broad ring geometry in 
an electromagnetic cavity with quantized photon modes 
remains unexplored beyond the Markovian approxima- 
tion. 

In the present work, we explore the transient effects 
of electronic transport through a broad quantum ring 
in a linearly polarized electromagnetic cavity coupled to 
electrically biased leads. This electron-photon coupled 
system under investigation can be manipulated by tun- 
ing the applied magnetic field and the polarization of the 
photon field. A time-convolutionless (TCL) version of 
the GME is utilized to project the time evolution onto 
the central system by taking trace with respect to the op- 
erators in the leads. 48 50 We demonstrate the transient 
transport properties by showing the many-body (MB) en- 
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ergy spectra, the time-dependence of the electric charge 
and current, the magnetic-field dependence of the total 
charge current with (w) or without (w/o) photon cavity, 
the charge density distribution, the normalized current 
density distribution and the local current coming from 
an occupation redistribution of the MB states in the cen- 
tral quantum ring system. 

The paper is organized as follows. In Sec. II, the the- 
oretical model is described. The electron system is em- 
bedded in an electromagnetic cavity by coupling a many- 
level electron system with photons using the full photon 
energy spectrum of a single cavity mode. In Sec. Ill, we 
show the numerical results for the dynamical transient 
transport properties for different magnetic field and pho- 
ton field polarization. Concluding remarks will be pre- 
sented in Sec. IV. 
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FIG. 1. (Color online) Equipotential lines in the central ring 
system connected to the left and right leads. Note that the 
isolines are refined close to the bottom of the ring structure. 



II. MODEL AND THEORY 

In this section, we describe the central system poten- 
tial Vs for the broad quantum ring and its connection to 
the leads. The electronic ring system is embedded in an 
electromagnetic cavity by coupling a many-level electron 
system with photons using the full photon energy spec- 
trum of a single cavity mode. The central ring system is 
described by a MB system Hamiltonian Hs with a uni- 
form perpendicular magnetic field, in which the electron- 
electron interaction and the electron-photon coupling to 
the x- or y-polarized photon field is explicitly taken into 
account. We employ the TCL-GME approach to explore 
the non-equilibrium electronic transport when the system 
is coupled to leads by a transient switching potential. 



A. Quantum ring connected to leads 

The system under investigation is a broad quantum 
ring connected to left and right leads / £ {L, R} with 
identical parabolic confining potentials 



TABLE I. Parameters of the central region ring potential. 



Vi(r) = -m*tf y 2 



(1) 



in which the characteristic energy of the confinement is 
hQo = 1.0 meV and m* = 0.067m e is the effective mass 
of an electron in GaAs-based material. 

The quantum ring is embedded in the central system of 
length L x = 300 nm situated between two contact areas 
that will be coupled to the external leads, as is depicted 
in Fig. 1. The system potential is described by 
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B. Central system Hamiltonian 

The time-evolution of the closed system with respect 
to t = 



U s (t) = exp[--H s t 



is governed by the MB system Hamiltonian 51 
I (h 



(3) 



d 2 r ffi(r) 



2m* 
+ H K 



V+- 

i c 



A(r) + A ph (r) 



(4) 



The first term includes a constant magnetic field B = 
Bz, in Landau gauge being represented by A(r) = 
—Byx. The second term is the exactly treated electron- 
electron interaction 



H ee = 
where 



d 2 r 



J d 2 r^t( r )^t (r / ) y ee(r5r / ) ^ (r / ) ^ (r ) 5 (5) 



Vs(*) = ^2 V i ex P - (Pxi(x ~ X )) 2 ~ (PyiVY 
i=l 

1 — 2 2 



-m*n o y 
with parameters from table I. 



V ee (r,r') 



_/|2 



(6) 



(2) with e > being the magnitude of the electron charge 
and rj = 1.0 x 10 -15 nm being a numerical regularization 
parameter. In addition, the last term in Eq. (4) indi- 
cates the quantized photon field, where a and a) are the 



3 



photon annihilation and creation operators, respectively, 
and Huj is the photon excitation energy. The photon field 
interacts with the electron system via the vector poten- 
tial 



A ph (r) = A(a + a f ) 



ex 

G7, 



TEqh 
TE101 



(7) 



for longitudinally-polarized (x-polarized) photon field 
(TEqh) or transversely-polarized (y-polarized) photon 
field (TE101). The electron-photon coupling constant 
gEM _ e Aa w Q w / c scales with the amplitude A of the 
electromagnetic field. For reasons of comparison, we also 
consider results without photons in the system. In this 
case, A ph (r) and huoa^a drop out from the MB system 
Hamiltonian in Eq. (4). 



C. Time-convolutionless generalized master 
equation approach 

The TCL-GME 50 is an alternative non-Markovian 
master equation to the Nakajima-Zwanzig (NZ) 
equation, 52 55 which is local in time. We assume, the 
initial total statistical density matrix can be written as a 
product of the system and leads density matrices, before 
switching on the coupling to the leads, 



W(0) =PL®PR®PS(0), 



(8) 



with pi, I E {L : R} : being the normalized density matri- 
ces of the leads. The coupling Hamiltonian between the 
central system and the leads reads 



Hr 



it) = E / ^ [&b)c q i + C\^{q) 

7 r d J 



l=L,R l 



(9) 



Here, C ql is the electron creation operator for state q and 
lead I and 

a/3 a 

with the creation operator, for the single-electron 
state (SES) a in the central system, i.e. the eigenstate a 
of the first term of Eq. (4) with A ph (r) = 0. The coupling 
is switched on at t = via the switching function 



x\t) = 1 



1 



(11) 



with switching parameter a 1 . Eq. (10) is written in the 
system Hamiltonian MB eigenbasis {|c*)}. The coupling 
tensor 51 



d 2 r / dV ^(r')5i g (r,r')Vf(r) (12) 



couples the extended lead SES {i/j q i(r)} with energy spec- 
trum {e l (q)} to the system SES {^f(r)} with energy 



spectrum {E a } that reach into the contact regions, 54 Q l s 
and 11/, of system and lead I, respectively, and 



g l aq (r, r>) = g l exp [S l x (x - x'f - 5 l y (y - y'f] 
\E a -e l (q)\' 



x exp 



(13) 



Here, g l is the lead coupling strength. In addition, S l x 
and 5y are the contact region parameters for lead I in x- 
and y-direction, respectively. Moreover, denotes the 
affinity constant between the central system SES energy 
levels {E a } and the lead energy levels {e l (q)}. 

In this work, we derive the TCL-GME 50 in the 
Schrodinger picture. In this picture, the reduced den- 
sity operator (RDO) of the system, 



p s (t) = Tr L Tr R \W(t)], 



(14) 



evolves to second order in the lead coupling strength in 
time via 



■j:[Hs,Ps{t)] 



l=L,R 



/ dq k l (q),Cl l (q,t)p s (t) 



with 



-f(e l (q)){ps(t),&(q,t)}] + H.c. 



n l (q,t) = ^ X l (t)e X p(-^te l (q) 
xU s (t)tl l (q,t)Ul(t), 



(15) 



tl l (q,t) = [ dt' . 
Jo 



exp ( -t'e l (q) ) X l (t>) 



(16) 



(17) 



and f(E) being the Fermi distribution function. 

Comparing this equation to the corresponding NZ 
equation, 52 55 

pT(t) = -l[H s ,pT(t)} 



J2 I dq[% 1 (q),ti l (q,t)}+ H.C. 
il=L,R 



(18) 



with 
Cl l (q,t) 



= ^X l Ws(t) J*dtf 
X X l (t')U l (q,t')]uUt) 



exp I -(t'-t)e l (q) 



(19) 



and 



tl l ( q ,t') = uUt')[± l Hq)pT(t') 

-/(^^{pTin^iq^Mn (20) 
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we note that we reobtain the TCL equation, if we set 



pr(t') = ui(t-t')p s (t)Us(t-t'), 



in Eq. (20) (which enters the kernel of Eq. (18)), but let 
P<P(*) = Ps(t) in the first term of Eq. (18). In other 
words, in the Schrodinger picture, the NZ kernel takes 
the central system time propagated RDO (which lets it 
become convoluted), while the TCL kernel takes just the 
unpropagated RDO. The deviation between the two ap- 
proaches is therefore only of relevance when the central 
system is far from a steady state and when the coupling to 
the leads is strong. It is our experience that the positivity 
conditions 56 for the MB state occupation probabilities in 
the RDO are satisfied to a higher system- lead coupling 
strength in the TCL case. The more involved quantum 
structure demands a stronger system-lead coupling than 
in our earlier work. 51 The numerical effort of the two ap- 
proaches is of similar magnitude. Both cases allow for 
a t-independent inner time integral over t\ which can 
be integrated successively with increasing t (increasing 
integration domain). 57 The RDO is inside (NZ) or out- 
side (TCL) of the inner time integral, but the required 
number of matrix multiplications is equal. 



III. NON-EQUILIBRIUM TRANSPORT 
PROPERTIES 

In this section, we investigate the non-equilibrium elec- 
tron transport properties through a quantum ring sys- 
tem, which is situated in a photon cavity and weakly 
coupled to leads. We assume GaAs-based material with 
electron effective mass m* = 0.067m e and background 
relative dielectric constant n = 12 A. We consider a 
single cavity mode with fixed photon excitation energy 
huo = 0.4 meV. The electron-photon coupling constant in 
the central system is g EM = 0.1 meV. Before switching 
on the coupling, we assume the central system to be in 
the pure initial state with electron occupation number 
^Ve,init = an d photon occupation number N p h,\mt — 1 
of the electromagnetic field. 

An external perpendicular uniform magnetic field is 
applied through the central ring system and the lead 
reservoirs. The area of the central ring system is A « 
2 x 10 4 nm 2 so that the magnetic field corresponding to 
the flux quantum <E> is B = $o/A « 0.2 T. The temper- 
ature of the reservoirs is assumed to be T = 0.5 K. The 
chemical potentials in the leads are hl = 2 meV and 
jiR = 0.9 meV leading to a source-drain bias window 
A/i = 1.1 meV. To facilitate inelastic scattering pro- 
cesses between the SES in the central ring system and 
the SES in the lead /, we allow for coupling of highly en- 
ergetically different states by letting the affinity constant 
A^ = 4.0 meV. 58 In addition, we let the contact region 
parameters for lead I £ {L, R} in x- and y-direction be 
5 l x = 5 l = 4.39 x 10 -4 nm -2 . The system-lead coupling 



There are several relevant length and time scales that 
should be mentioned. The two-dimensional magnetic 



(21) length is I = [ch/(eB)\ 



1/2 



25.67[£(T)]- 1 / 2 nm. The 



ring system is parabolically confined in the y-direction 
with characteristic energy Hflo = 1-0 meV leading to a 
modified magnetic length scale 



ft 



1/2 



ra*^o/ ^1 + {eB/(m*cQ )} 2 
33.74 



^/1 + 2.982[£(T)] 2 



nm. 



(22) 



Correspondingly, the system-lead coupling strength is 
then g l dw 2 = 39.85 meV for magnetic field B = 0.1 T 
and g l al /2 = 38.22 meV for magnetic field B = 0.225 T. 
The time-scale for the switching on of the system-lead 
coupling is (a 1 ) -1 — 3.291 ps, the single-electron state 
(1ES) charging time-scale ties ~ 30 ps, and the two- 
electron state (2ES) charging time-scale T2es ^ 200 ps 
described in the sequential tunneling regime. We study 
the transport properties for < t < T2es, when the sys- 
tem has not yet reached a steady state. 

In order to understand the non-equilibrium dynamical 
behavior of the charge distribution in the system, we de- 
fine the time-dependent magnitude of charge on the left 
part (x < 0) of the ring 



} L s (t)= I dx I dyp(v,t) (23) 



and the time-dependent magnitude of charge on the right 
part (x > 0) of the ring 

Q§(t)= / dx / dyp(r,t). (24) 

The space- and time-dependent charge density, 

p(r,t) = Tr\ps(t)p(r)], (25) 

is the expectation value of the charge density operator 

p(r) = e^ f (r)Vi(r). (26) 

In order to explore the magnetic field influence on the 
charge currents from and into the leads, we define the 
charge current from the left lead into the system by 



I L (t) = Tv[p L s (t)Q] 



(27) 



strength g l = 0.2058 meV/nm 



3/2 



Here, Q = eN is the charge operator with number oper- 
ator TV and the time-derivative of the RDO in the MB 
basis due to the coupling to the lead I £ {L, R} 



-H.c (28) 



Ps(t) = Jdq [% l (q),[n l (q,t)ps(t) 
f(6 l (q)){p s (t),n l ( q ,t)} 
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Similarly, the charge current from the system into the 
right lead can be expressed as 



I R (t) = -Tv[^(t)N}. 



(29) 



To get more insight into the local current flow in the ring 
system, we define the top local charge current through 
the upper arm (y > 0) of the ring 



itopW 



dy j x (x = 0,y,t) 



(30) 



and the bottom local charge current through the lower 
arm (y < 0) of the ring 



^bottom (t) = / dy j x (x = 0,y,t) . 

J — OO 



Here, the charge current density, 



(31) 



(32) 



is given by the expectation value of the charge current 
density operator, 



j(r) = j p (r) +j d (r), 



(33) 



decomposed into the paramagnetic charge current den- 
sity operator, 



eh 



2mi 



(34) 



^(r)(V^(r)) - (V^(r))^(r) 

and the diamagnetic charge current density operator, 

Ur)=}THr)+}f(r). (35) 
The latter consists of a magnetic component, 



jmag (r) = f_ A(r) ^t (r) ^ (r)> 



rn 



and photonic component, 



jP\r) = -AP h (r)^ (r )^( r ). 

rn 



(36) 



(37) 



Furthermore, to understand better the driving schemes 
of the dynamical transport features, we define the total 
local charge current 



hl(t) = itop(^) + ^bottom (t) 

and circular local charge current 

Icl(t) = ^(hop(t) _ ^bottom (t)). 



(38) 



(39) 



Below, we shall explore the influence of the applied mag- 
netic field and the photon fiel polarization on the non- 
equilibrium quantum transport in terms of the above 
time-dependent charges and currents in the broad quan- 
tum ring system connected to leads. 
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FIG. 2. (Color online) MB energy spectrum of system Hamil- 
tonian Hs versus magnetic field B in units of tesla (T). The 
states are differentiated according to their electron content 
iV e : zero-electron states (iV e = 0, OES, green dots), single 
electron states (A^ e = 1, 1ES, red dots) and two electron states 
(N e = 2, 2ES, blue crosses). The photon field is x-polarized. 



A. Photons with x-polarization 

In this subsection, we focus on our results for x- 
polarized photon field. Figure 2 shows the MB energy 
spectrum of the system Hamiltonian Hs including the 
electron-electron and electron-photon interactions. The 
MB-energy levels are assigned different colors according 
to their electron content N e . For N e = (green dots), 
the MB states differ in energy by multiples of the pho- 
ton energy huj according to their photon content N p h 
independently of the applied magnetic field. The bias 
window contains a number of SES (red dots) of which 
the two lowest ones have the highest occupation. The 
SESs show crossing behavior at half integer flux quanta. 
The state- dependent effective magnetic flux ring area al- 
lows for small variations of the crossing period Bq. The 
crossings at integer flux quanta are usually avoided due 
to the ring rotation symmetry violation, which is mainly 
coming from the presence of the central system contact 
regions to the leads. In general, the SES and, in partic- 
ular, the two electron states (2ES, blue crosses) tend to 
increase in energy with higher magnetic field. 

Figure 3 illustrates the central region charging of 1ES 
and 2ES as a function of time with the initial condi- 
tions Qies(O) = Q2Es(0) = since we selected an initial 
state with AT e?init = 0. In the low magnetic field regime, 
in the case B = 10 -5 T, we notice that (t, Qies) = 
(200ps,0.809e) and (t, Q 2 es) = (200 ps, 0.055e). In the 
high magnetic field regime, in the case B = 1 T, we 
notice that (t, Qies) = (200 ps, 0.911e) and (t, Q2ES) = 
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FIG. 3. (Color online) Switching function x l (t) (solid red), 
charge of all 1ES for B = 10" 5 T (dashed green) and 
B = 1.0 T (dotted blue), and charge of all 2ES for B = 10" 5 T 
(dotted purple) and B = 1.0 T (dash-dotted cyan) as a func- 
tion of time. The photon field is x-polarized. 

(200 ps, 0.012e). In general, the 2ES are occupied slower 
than the 1ES indicating the sequential tunneling pro- 
cesses. But more importantly, the energetic location of 
the 2ES and their shifting above the bias window by the 
Coulomb interaction attenuates the 2ES occupation. As 
is shown in Fig. 2, the magnetic field plays a role to 
increase further the energy difference of the 2ES with 
respect to the 1ES, thus enhancing the 1ES occupation 
by £Qies = 0.102e while reducing the 2ES occupation 
by SQ2ES — — 0.043e. We note that the earlier men- 
tioned time-scales (a 1 ) -1 = 3.291 ps, ties ~ 30 ps and 
t 2es ^ 200 ps are in agreement with Fig. 3. 
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FIG. 4. (Color online) The left charge current II (solid red) 
and the right charge current Ir (long-dashed green) versus 
the magnetic field with (w) x-polarized photon field at t — 
200ps. For comparison: left charge current II (short-dashed 
blue) and right charge current Ir (dotted purple) in a purely 
electronic central system, i.e. without (w/o) photon cavity. 



In Fig. 4, we show the current from the left lead into 
the ring system II (solid red curve) and the current from 
the ring system to the right lead Ir (long-dashed green 
curve) as a function of magnetic field at time t = 200 ps. 
The similar values of Il{B) and Ir(B) indicate that 
the short-time regime charging of 1ES is nearly com- 
pleted at t = 200 ps meaning that the total charging 
has slowed down by more than an order of magnitude. 
Moreover, we see clear oscillations of the current with pe- 
riod Bq « 0.2 T: the first minimum current at B = 0.1 T 
corresponds to the situation of a half flux quantum, in 
which the left charge current II = 0.273 nA and the right 
charge current Ir = 0.261 nA and the maximum current 
at B = 0.225 T is corresponding to the case of one flux 
quantum, in which the left charge current II = 1.183 nA 
and the right charge current Ir = 1.168 nA. These ob- 
servations are in agreement with the Aharonov-Bohm 
(AB) oscillations of the steady state 12 14 but superim- 
posed and modified by electron-electron correlation ef- 
fects and the non-equilibrium situation. In addition, the 
electron-photon coupling suppresses the constructive in- 
terference of AB-phases in the integer flux quantum sit- 
uation as can be seen from a comparison with the purely 
electronic system results in Fig. 4 (short-dashed blue and 
dotted purple curve). 

Figure 5 shows the time-evolution of the left total cur- 
rent the right total current the top local 
current 7 top (t), and the bottom local current /bottom(£)- 
In the case of B = 0.1 T (a half flux quantum) shown 
in Fig. 5(a), the maximum value of the charge current 
from the left lead at t = 15.30 ps is I L = 4.190 nA. In 
addition, the minimum value of the charge current to the 
right lead at t = 12.25 ps is Ir = —0.511 nA. The nega- 
tive right charge current indicates that the central system 
is also charged from the right lead during the transient 
phase. In the case of B = 0.225 T (one flux quantum) 
shown in Fig. 5(b), the maximum value of the current 
from the left lead at t = 15.35 ps is II = 4.334 nA; 
and the minimum right charge current at t = 10.80 ps is 
Ir = —0.387 nA. It is then interesting to realize that the 
magnetic field enhances the charge accumulation from 
the left lead, while it suppresses the short-time regime 
charging from the right lead. Hence, the integer mag- 
netic flux quantum case assists the net current flow al- 
ready in the highly non-equilibrium situation in the very 
beginning. The local charge transport may differ in di- 
rection in the ring arms due to the persistent magnetic 
field induced ring current. 

In Fig. 6, we illustrate the normalized charge current 
density vector field j(r,t) in the central quantum ring 
system in the case of x-polarized photon field in the long- 
time response regime t = 200 ps, i.e. when the 2ES get 
charged. For magnetic field B = 0.1 T, a clear counter- 
clockwise vortex located close to the left lead can be 
found dominating the current flow pattern in the central 
ring system as shown in Fig. 6(a). 16 The vortex circu- 
lation direction is in agreement with the Lorentz force, 
since the enclosed area is threaded by much less than half 
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FIG. 5. (Color online) Local current through the top ring 
branch It op (t), local current through the bottom ring branch 
/bottom (t), current from the left lead into the system iz,(t), 
and current from the system into the right lead /#(£) as a 
function of time for (a) B = 0.1 T and (b) B = 0.225 T in 
the case of x-polarized photon field. 



a flux quantum. Due to the geometrical position of the 
vortex and the current continuity condition, clockwise 
current direction is favored for the ring system. However, 
the counter-clockwise vortex appears relatively weak for 
magnetic field B = 0.225 T present at both left and right 
lead connection area as shown in Fig. 6(b), while the to- 
tal local current through the whole central system from 
the left to the right lead is large. Additionally, for a 
later comparison with the y-polarized photon field, Fig. 
6(c) shows the current density for B = 0.425 T (two flux 
quanta), which is similar to Fig. 6(b) (one flux quanta) 
with the vortex circulation on both sides being slightly 
more significant. 

Figure 7 illustrates the time-dependent charge on the 
left part of the ring Q$ (t) and the time-dependent charge 
on the right part of the ring Q§(t). In the case of 
B = 0.1 T shown in Fig. 7(a), both and Q§ are 
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FIG. 6. (Color online) Normalized charge current density 
vector field in the central system for (a) B = 0.1 T, (b) B = 
0.225 T, and (c) B = 0.425 T at t = 200 ps in the case of 
x-polarized photon field. 



increasing almost monotonically in time. In the long- 
time response regime t = 200 ps, Qg(t) = 0.742e is much 
higher than Q§(i) = 0.234e. This implies charge accu- 
mulation mainly on the left hand side of the quantum 
ring in the case of magnetic field with half integer flux 
quantum with enhancement of the electron dwell time 
on the left-hand side of the ring and suppression of the 
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FIG. 7. (Color online) Charge in the left (Qs(t)) or right half 
(Qs (£)) of the central quantum ring system as a function of 
time for (a) 5 = 0.1 T and (b) B = 0.225 T. The photon field 
is x-polarized. 

electron dwell time on the right-hand side. 

In the case B = 0.225 T shown in Fig. 7(b), both Q£ 
and Q§ exhibit oscillations in time after the short-time 
charging regime. This implies that the charge accumula- 
tion manifests itself in oscillating behavior between the 
left and the right part of the quantum ring in the case of 
magnetic field with integer flux quantum. The oscillation 
amplitude is decreasing in time due to the dissipation ef- 
fects caused by the coupling to the leads. In the long-time 
response regime t = 200 ps, Q$(t) = 0.423e is of similar 
magnitude than Q§(t) = 0.446e, which is by difference 
to the half integer flux quantum case. It is interesting to 
notice that the oscillating period of the charges is around 
r = 100 ps corresponding to a characteristic energy 
5Eq « 0.04 meV. The MB energies of the mostly occu- 
pied MB levels are Ef = 1.4038 meV E§ = 1.3664 meV 
such that AEg 10 = 0.0374 meV. The corresponding two- 
level (TL) oscillation period of the closed system would 
be 

t tl = HI ps. In the non-equilibrium open system, 



the TL oscillation period is r^ L = 94 ps or r R L = 100 ps 
when we take the time intervals between the first and 
second maxima of Q$(t) and Q§(t), respectively. The 
full numerical calculation including all MB levels shown 
in Fig. 7(b) yields the left and right charge oscillation 
period, r L = 96 ps and r R = 110 ps, respectively. The 
system is far from equilibrium at the earlier maximum, 
thus reducing in particular the left period t^ l with re- 
spect to r^ L . However, we find that also the other MB 
states change the periods when comparing r L with t^ l 
and r R with r^ L . 

We would like to bring attention to the fact that charge 
balances like Ql = II — h\ an d Qr = hi — Ir would not 
be satisfied. This is because the SES that are filled from 
the left lead or emptied to the right lead are in general 
not restricted to a single half of the central system, but 
extended over the whole system. 

Figure 8 shows the charge density distribution in the 
central quantum ring system in the case of x-polarized 
photon field with the magnetic field (a) B = 0.1 T, (b) 
B = 0.225 T, and (c) B = 0.425 T at t = 200 ps. In 
the case of B = 0.1 T (half flux quantum) shown in Fig. 
7(a), the electrons are highly accumulated on the left- 
hand side of the quantum ring with very weak coupling 
to the right lead, and hence strongly blocking the left 
charge current and suppressing the right charge current, 
as it was shown previously in Fig. 4 (marked by the up- 
arrow). For half integer flux quantum, the electron dwell 
time on the left-hand side of the ring is enhanced relative 
to the electron dwell time on the right-hand side of the 
ring due to destructive phase interference on the right 
hand side. The density accumulates then mainly on the 
left hand side of the ring forming a long-living state and 
the magnetic field evoked vortex on the right hand side 
contact area is suppressed. 

In the B = 0.225 T case (one flux quantum) shown in 
Fig. 8(b), the electrons manifest oscillating feature be- 
tween the left and right end of the quantum ring. At 
time t = 200 ps, the electrons are nearly equally well 
accumulated on both sides of the quantum ring. This 
phenomenon is related to the manifestation of current 
peaks observed in Fig. 4 (marked by the down-arrow). 
The charge distribution is rearranged when compared to 
magnetic field B = 0.1 T. We observe a depletion of 
about 50% at the left-hand contact region with equivalent 
charge augmentation on the right-hand contact region. 
The magnetic field B = 0.225 T with integer flux quan- 
tum enhances the likelihood for electrons to flow through 
the quantum ring to the right-hand side of the central 
system and further to the right lead. Additionally, Fig. 
8(c) shows the charge density for B = 0.425 T (two flux 
quanta), which is similar to Fig. 8(b) (one flux quantum). 

In Fig. 9, we show the magnetic field dependence of 
the partial local currents I top and /bottom through the 
top and bottom arms, the total local current I t \ across 
x = 0, and the circular local current J c i, which are con- 
venient tools to study the relative importance of local 
"persistent" current flows induced by the magnetic field 
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FIG. 8. (Color online) Charge density distribution 
p(r, t) (e/a^) in the central system for (a) B — 0.1 T, (b) 
B = 0.225 T, and (c) B = 0.425 T in the x-polarized photon 
field case at t = 200 ps. 



in the long-time response transient time regime. We aver- 
aged the local currents over the time interval [180, 220] ps 
around t = 200 ps to soften possible high frequency fluc- 
tuations (compare with Fig. 5). In general, the top local 
current exhibits opposite sign to the bottom local cur- 
rent, and hence the circular local current (solid purple) 
is usually larger in magnitude than the total local cur- 
rent (dotted blue). The local current through the two 
current arms, I t \, is strongly suppressed in the case of 
half integer flux quanta showing a very similar behav- 
ior to the nonlocal currents II and Ir (Fig. 4). This is 
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FIG. 9. (Color online) Local current through the top arm 
of the ring I top (solid red), local current through the bottom 
arm of the ring /bottom , total local current It\ , and the circular 
local current I c \) versus the magnetic field averaged over the 
time interval [180, 220] ps in the case of x-polarized photon 
field. 



because the destructive interference in the quantum ring 
enhances the back scattering for magnetic flux with half 
integer quanta. 

In the absence of magnetic field 5 = 0, the circular 
current is identical to zero due to the symmetric situa- 
tion for both ring arms. Moreover, it is interesting to 
note that the circular local current I c \ reaches 1.347 nA 
for less than half a flux quantum (at B = 0.05 T), in- 
creases further until B = 0.45 T with a maximum value 
max |/ c i| = 2.844 nA and decreases again for B > 0.45 T. 
The magnetic component of the diamagnetic part of the 
circular local current increases linearly with the magnetic 
field 5, but the paramagnetic part guarantees a behavior, 
which is closer to being periodic with the flux quantum. 
The periodic part is in analogy to the behavior for a ring 
of infinitesimal width. 59 In the case of high magnetic field 
regime (B > 0.45 T), a comparison with Fig. 2 shows 
that the different flux periods of different MB-states in 
the finite-width ring lead to destructive interference ef- 
fects reducing the periodic oscillations considerably. The 
most common direction of the circular current is clock- 
wise. This is because the vortices at the lead connection 
areas are threaded by less than half a flux quantum lead- 
ing to counter-clockwise vortex rotation direction. Then, 
as a consequence of charge continuity and the vortex lo- 
cation outside the ring radius, clockwise direction for I c \ 
is preferred. 



B. Photons with y-polarization 

In this subsection, we focus on the y-polarized pho- 
ton field situation and compare with the results for the 
x-polarized photon field. Figure 10 shows the MB en- 
ergy spectra of the system Hamiltonian Hs in the case 
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FIG. 11. (Color online) Switching function x\t) (solid red), 
charge of all 1ES for B = 10" 5 T (dashed green) and 
B = 1.0 T (dotted blue), and charge of all 2ES for B = 10" 5 T 
(dotted purple) and B = 1.0 T (dash-dotted cyan) as a func- 
tion of time. The photon field is y-polarized. 



FIG. 10. (Color online) MB energy spectrum of the sys- 
tem Hamiltonian Hs versus magnetic field B within the bias 
window energy range for (a) x-polarized and (b) y-polarized 
photon field. The states are differentiated according to their 
electron content N e : zero-electron states (N e = 0, OES, green 
dots) and single electron states (N e = 1, 1ES, red dots). 

of (a) x-polarized and (b) y-polarized photon field. We 
note in passing that Fig. 10(a) magnifies a part of the 
MB spectrum of Fig. 2. The mostly occupied levels are 
the two levels around 1.4 meV. In the cases of both x- 
and y-polarized photon field, we see the MB energy de- 
generacy around B = 0.1 and 0.325 T related to the de- 
structive AB phase interference. However, in the case of 
y-polarization, an extra MB energy degeneracy is found 
at B — 0.425 T. This degeneracy is related to a photon 
suppressed current dip, i.e. not related to AB oscilla- 
tions. 

Figure 11 illustrates the charge of 1ES and 2ES as a 
function of time. In the case of low magnetic field regime 
B = 10" 5 T, we notice that (t, Qies) = (200 ps, 0.847e) 
and (£, Q2ES) = (200 ps, 0.055e). In the case of high mag- 
netic field regime B = 1 T, we notice that (£,Qies) = 
(200ps,0.891e) and (t, Q 2 es) = (200 ps, 0.017e). The 
2ES are occupied much slower than the 1ES due to 
similar reasons than for the x-polarized photon field. 
The magnetic field enhances the 1ES occupation by 
^Qies = 0.044e while it suppress the 2ES occupation 
by £Q2ES = — 0.038e. Therefore, in comparison with the 
case of x-polarization shown in Fig. 3, we realize that 
the y-polarized photon field mildens the influence of the 
applied magnetic field on the charging feature. 

Figure 12 shows the left charge current II (solid red) 
and the right charge current Ir (dashed green) as a func- 
tion of magnetic field in the case of y-polarized photon 
field at t = 200ps. The similar values of and Ir(B) 
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FIG. 12. (Color online) Left charge current II (solid red) and 
right charge current Ir (long-dashed green) versus the mag- 
netic field with (w) y-polarized photon field at t = 200ps. For 
comparison: left charge current II (short-dashed blue) and 
right charge current Ir (dotted purple) in a purely electronic 
central system, i.e. without (w/o) photon cavity. 

agree well with the long-time response regime slow-down 
in charging predicted in Fig. 11, which is almost com- 
pleted for the 1ES. Charge current oscillations are shifted 
slightly from the period Bo w 0.2 T due to the broad 
ring geometry. Moreover, the oscillation amplitude and 
extrema positions show more unexpected features than 
in the case of x-polarized photon field. The first current 
minimum is at magnetic field B = 0.1 T (with a half 
flux quantum) with left charge current II = 0.303 nA 
and right charge current Ir = 0.298 nA. At magnetic 
field B = 0.225 T corresponding to the case of one flux 
quantum, the left charge current II = 1.007 nA and 
the right charge current Ir = 0.935 nA. It is interest- 
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ing to point out that the magnetic field dependence of 
the charge current exhibits a pronounced dip at mag- 
netic field B = 0.425 T (two flux quanta) in the case of 
y-polarized photon field that is not present in the case of 
x-polarized photon field. 

The dip structure in the charge current at B = 0.425 T 
is due to the above mentioned degeneracy of the MB 
energy spectrum, which strongly suppresses the photon- 
assisted tunneling feature. Furthermore, the charge cur- 
rent can be enhanced by the y-polarized photon field at 
magnetic field with half integer flux quantum: the y- 
polarized photon field significantly influences the quan- 
tum interference of the circular local current flow includ- 
ing the destructive interference feature of the charge cur- 
rent. 
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FIG. 13. (Color online) Local current through the top ring 
branch (itop(£)) and bottom ring branch (/bottom (£)), current 
from the left lead (Ic(t)) and into the right lead (/#(£)) for 
(a) B = 0.1 T and (b) B = 0.225 T in the case of y-polarized 
photon field. 

Figure 13 illustrates the time-evolution of the left total 
current Jz,(t), the right total current Jr(£), the top local 
current 7 top (£), and the bottom local current /bottom (t). 



In the case of B = 0.1 T shown in Fig. 13(a), the max- 
imum value of the current from the left lead into the 
system at t = 15.55 ps is /l(£) = 4.303 nA. Furthermore, 
the minimum value of the charge current into the right 
lead at t = 12.00 ps is I R = -0.494 nA. The negative 
right charge current indicates that the central system 
is charged from the left and the right lead for a short 
time. This charging from the right is a little weaker than 
in the case of x-polarized photon field. In the case of 
B = 0.225 T shown in Fig. 5(b), the maximum value of 
the left current at t = 15.54 ps is II = 4.446 nA and 
the minimum value of the right current at t = 10.80 ps 
is Ir = —0.337 nA. Hence, the integer magnetic flux en- 
hances the charge accumulation from the left lead, while 
it suppresses the short-time regime charging from the 
right lead assisting the net current flow from the left to 
the right already in the highly non-equilibrium situation 
in the very beginning. The local charge transport may 
differ in direction in the ring arms due to the "persistent" 
current induced by the magnetic field. 

In Fig. 14, we demonstrate the normalized charge cur- 
rent density vector field j(r, t) in the central ring system 
for the magnetic field, (a) B = 0.1 T, (b) B = 0.225 T, 
and (c) B = 0.425 T, in the long-time response regime 
t = 200 ps in the case of y-polarized photon field. For 
magnetic field B = 0.1 T, a clear counter-clockwise vor- 
tex can be found being associated with a long-living lo- 
calized state which is strongly dominating the current 
flow pattern in the central ring system, as is shown in 
Fig. 14(a). However, for magnetic field B = 0.225 T, 
this counter-clockwise vortex appears weaker relative to 
the total local current, but is present at both contact re- 
gions as shown in Fig. 14(b). Figures 14(a) and 14(b) 
are similar to Figs. 6(a) and 6(b) meaning that the local 
current flow is mainly governed by AB interference with 
the photon polarization having only a minor effect. 

Figure 14(c) shows the current density for B = 0.425 T 
(two flux quanta), which is similar to Fig. 14(a) (a half 
flux quantum) and not to the one flux quantum case as 
for x-polarization (similarity between Fig. 6(c) and Fig. 
6(b)). This feature is not predicted by the AB effect, but 
is caused by the influence of the y-polarized photons. 
However, the impact of a MB spectrum degeneracy of 
the mostly occupied MB states (Fig. 10(b)) on the local 
current flow structure is similar whether the degeneracy 
is in agreement with the AB effect (Fig. 14(a)) or not, 
i.e. originates from the photons (Fig. 14(c)). 

Figure 15 shows the time-evolution of Qg(i) and 
Q§(i). In the short-time response regime at t — 65 ps, 
the charge on the left and right part of the ring are 
Q L S = 0.548e and Q§ = 0.318e for B = 0.225 T, re- 
spectively; the charge on the left and right part of the 
ring are = 0.610e and Q§ = 0.279e for B = 0.1 T. 
Consequently, the electron dwell time on the left-hand 
side of the ring is enhanced while the electron dwell time 
on the right-hand side is suppressed for y-polarized pho- 
ton field and integer flux quantum; this feature is though 
more pronounced for a half flux quantum and y-polarized 
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FIG. 14. (Color online) Normalized charge current density 
vector field in the central system for (a) B = 0.1 T, (b) B = 
0.225 T and (c) B = 0.425 T at t = 200 ps in the case of 
y-polarized photon field. 



photon field, but was already described for x-polarized 
photons (Fig. 7(a)). The reason for the difference in left 
and right dwell time in the integer flux quantum case is 
a low frequency oscillation of the most important close- 
in-energy levels. In the long-time response regime at 
time t = 200 ps, the picture is very similar to the x- 
polarized photon field case: for B = 0.225 T, the left 




FIG. 15. (Color online) Charge in the left half (Qs(t)) or 
right half (Qs (£)) of the central system as a function of time 
for (a) B = 0.1 T and (b) B = 0.225 T. The photon field is 
y-polarized. 



and right charges are of similar magnitude, Q$ = 0.462e 
and Q§ = 0.431e; and for B = 0.1 T, the charge is 
mainly accumulated at the left hand side, Q$ = 0.720e 
and Q§ = 0.256e. 

In the B = 0.225 T case, the MB energies of the 
mostly occupied MB levels are Ef = 1.3846 meV and 
Eg = 1.3683 meV such that A£^ 10 = 0.0163 meV. The 
energy level difference of the mostly occupied MB lev- 
els is only 44% of the case of x-polarized photon field: 
Eq 10 ~ 0.44 x Eg 10 . The corresponding TL oscillation 



period of the closed system would be r^ L = 254 ps. The 
oscillation period is too long to be observed clearly in Fig. 
15(b), but the first maximum of Qs(t) at t — 65 ps repre- 
sents the starting point of the low frequency oscillation, 
which is better visible in the TL system defined by the 
two mostly occupied states. Our findings suggest that the 
energy difference of the two mostly occupied levels con- 
trols not only the charge distribution, but also photonic 
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suppressions of the AB current. The different connec- 
tivity (probability density on the left or right ring part) 
to the leads found within the TL dynamic suggests that 
the probability of a photon coupled electron transition 
between these levels plays a major role in understanding 
the photonic modifications of the AB current pattern. 
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FIG. 16. (Color online) Charge density distribution 
p(r, t) (e/cfiu) in the central system for (a) B — 0.1 T, (b) 
B = 0.225 T, and (c) B = 0.425 T in the y-polarized photon 
field case at t = 200 ps. 

Figure 16 shows the charge density distribution in the 
central ring system in the case of y-polarized photon field 
for magnetic field (a) B = 0.1 T, (b) B = 0.225 T, and 
(c) B = 0.425 T at t = 200 ps. In the case of B = 0.1 T 
shown in Fig. 16(a), the electrons are highly accumulated 



on the left-hand side of the quantum ring with very weak 
coupling to the right lead, and hence strongly blocking 
the left charge current and suppressing the right charge 
current. 

In the case of B = 0.225 T shown in Fig. 16(b), the 
electrons manifest oscillating feature between the left and 
right end of the quantum ring. At time t = 200 ps, the 
electrons are equally well accumulated on both sides of 
the quantum ring. This situation is related to the man- 
ifestation of the current peaks observed in Fig. 12. The 
charge is redistributed when compared to magnetic field 
B = 0.1 T. We observe a depletion of about 50% on the 
left-hand side with equivalent charge augmentation on 
the right hand side. The magnetic field B = 0.225 T 
with integer flux quanta enhances the likelihood for elec- 
trons to flow through the quantum ring to the right-hand 
side of the central system and further to the right lead. 

Figure 16(c) shows the charge density for B = 0.425 T 
(two flux quanta), which is similar to Fig. 16(a) (a half 
flux quantum) and not to the one flux quantum case as 
for x-polarization (similarity between Fig. 8(c) and Fig. 
8(b)). This feature is not predicted by the AB effect, 
but is caused by the influence of the y-polarized pho- 
tons. However, the impact of a MB spectrum degener- 
acy of the mostly occupied MB states (Fig. 10(b)) on the 
density distribution is similar whether the degeneracy is 
in agreement with the AB effect (Fig. 16(a)) or not, i.e. 
originates from the photons (Fig. 16(c)). 
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FIG. 17. (Color online) Local current through the top ring 
arm (/top) and bottom ring arm (/bottom) and total local cur- 
rent (/ti) and circular local current (/ c i) versus the magnetic 
field and averaged over the time interval [180, 220] ps in the 
case of y-polarized photon field. 

In Fig. 17, we show the magnetic field dependence 
of the local currents / top and /bottom through the top 
and bottom arms, respectively, the total local current / t i 
across x = 0, and the circular local current I c \. We aver- 
aged the local currents over the time interval [180, 220] ps 
around t = 200 ps to soften possible high frequency fluc- 
tuations (compare with Fig. 13). In most cases (however 
less often than for x-polarized photons), the top local cur- 
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rent exhibits opposite sign to the bottom local current, 
and hence the circular local current (solid purple) is usu- 
ally larger than the total local current (dotted blue). The 
local current through the two current arms, J t i 5 is sup- 
pressed in the case of half integer flux quanta showing a 
similar behavior to the nonlocal currents II and Ir (Fig. 
12), but with more irregularities due to the stronger ef- 
fective influence of the y-polarized photon field. It is 
interesting to note that the current suppression dip at 
B = 0.425 T (marked by the blue arrow in Fig. 12) ap- 
pears also in the local current (blue dashed curve) flowing 
through both ring arms from the left to the right. 

The circular local current reaches a maximum absolute 
value of max 1 7ci| = 1.905 nA at B = 0.625 T, which is 
by 0.939 nA smaller than for x-polarization. It is clearly 
visible from a comparison of Fig. 17 and Fig. 9 that the 
circular current is considerably smaller than in the x- 
polarized photon case, while the total local current is of 
the same order. Thus, the capability of the magnetic 
field to drive a rotational ring current is weakened by 
having the electromagnetic field y-polarized. In particu- 
lar, this can be said about the diamagnetic part of the 
circular local current leading to the much smaller value 
Id = 0.675 nA at low magnetic field B = 0.05 T. The 
periodicity of the circular local current is preserved bet- 
ter for x-polarized photon field as is for the total local 
current. We note here in passing that it is not possible 
to understand the non-trivial magnetic field dependence 
of I c \ by resorting solely to a TL description. 

IV. CONCLUDING REMARKS 

We have presented a time-convolutionless generalized 
master equation formalism that allows us to calculate the 
non-equilibrium transport of Coulomb interacting elec- 
trons through a broad quantum ring in a photon cavity 
under the influence of a uniform perpendicular magnetic 
field. The topologically nontrivial broad ring geometry 
allows for substantial electron-electron correlations rela- 
tive to their kinetic energy and, hence, a large basis is 
required for sufficient numerical accuracy. The magnetic 
field, however, increases slightly the energy difference of 
the 2ES with respect to the 1ES, thus enhancing the 
1ES occupation while suppressing the 2ES occupation. 
The central quantum ring 1ES are charged quickly from 
both leads. Electron-electron correlation and sequential 
tunneling slow down the 2ES charging in the long-time 
response regime. Aharonov-Bohm charge current oscilla- 
tions can be recognized in the long-time response regime 
with magnetic field period Bo = ^o/A which is related 
to the flux quantum $ and ring area A. 

In the case of x-polarized photon field, we have found 
charge oscillations between the left and right part of the 
quantum ring when the magnetic field is associated with 
integer flux quanta. The oscillation frequency agrees 
well with the energy difference of the two mostly occu- 
pied states. The relatively high energy difference for x- 



polarized photons is related to a relatively high transient 
current through the ring. The amplitude of the charge 
oscillations through the quantum ring is decreasing in 
time due to dissipation effects caused by the coupling to 
the leads. In general, the local current through the up- 
per ring arm exhibits opposite sign to the local current 
through the lower ring arm. Hence, the "persistent" cir- 
cular local current is usually larger than the total local 
charge current through both ring arms from the left to 
the right. The persistent current shows a periodic behav- 
ior with magnetic field, but with a tendency to clockwise 
rotation due to the contact region vortex structure. 

In the case of y-polarized photon field, the magnetic 
field dependence of the left and right charge current ex- 
hibits a pronounced dip at magnetic field B = 0.425 T 
corresponding to two flux quanta that is therefore clearly 
not related to the Aharonov-Bohm effect. The dip is as- 
sociated with a degeneracy of the two mostly occupied 
1ES at magnetic field associated with two flux quanta. 
The additional level crossing appears only for y-polarized 
photons. The generally lower energy difference of the two 
mostly occupied MB states in the case of y-polarization 
disturbs the constructive phase interference condition for 
the bias driven charge flow through the quantum device 
and decreases the persistent current magnitude. 

In conclusion, we have demonstrated for our ring geom- 
etry that y-polarized photons disturb our system stronger 
than x-polarized photons, suppressing magnetic field in- 
duced currents and perturbing flux periodicity beyond 
finite width effects by enhancing or suppressing bias- 
driven currents. It is interesting to compare these find- 
ings to the quantum wire case, where it was found that 
mainly x-polarized photons attenuate the central system 
charging due to a closer agreement of the photon mode 
energy and the characteristic electronic excitation en- 
ergy in x-direction. 51 In this paper, we have considered 
a more complex geometry, which reduces effectively the 
y-confinement energy MIq = 1.0 meV. The characteristic 
electronic excitation energy in y-direction may therefore 
be much closer to the photon mode energy huo = 0.4 meV, 
thus leading to a relatively strong influence of the y- 
polarized photon field on the electronic transport. The 
conceived magnetic field influenced quantum ring system 
in a photon cavity could serve as an elementary quantum 
device for optoelectronic applications and quantum infor- 
mation processing with unique characteristics by control- 
ling the applied magnetic field and the polarization of the 
photon field. 
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